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ABSTRACT 

Context. The Swift satellite has given continuous data in the range 0.3-150 keV from s to 10 6 s for GRB060218 associated with 
SN2006aj. This GRB which has an unusually long duration (Tgo ~ 2100 s) fulfills the Amati relation. These data offer the opportunity 
to probe theoretical models for Gamma-Ray Bursts (GRBs) connected with Supernovae (SNe). 

Aims. We plan to fit the complete y- and X-ray light curves of this long duration GRB, including the prompt emission, in order to 
clarify the nature of the progenitors and the astrophysical scenario of the class of GRBs associated to SNe Ib/c. 
Methods. We apply our "fireshell" model based on the formation of a black hole, giving the relevant references. It is characterized 
by the precise equations of motion and equitemporal surfaces and by the role of thermal emission. 

Results. The initial total energy of the electron-positron plasma E'°J = 2.32 x 10 50 erg has a particularly low value similarly to the 
other GRBs associated with SNe. For the first time we observe a baryon loading B = 10~ 2 which coincides with the upper limit for 
the dynamical stability of the fireshell. The effective CircumBurst Medium (CBM) density shows a radial dependence n c t m oc r~ a with 
1.0 ^ a :£ 1.7 and monotonically decreases from 1 to 10~ 6 particles/cm 3 . Such a behavior is interpreted as due to a fragmentation in the 
fireshell. Analogies with the fragmented density and filling factor characterizing Novae are outlined. The fit presented is particularly 
significant in view of the complete data set available for GRB060218 and of the fact that it fulfills the Amati relation. 
Conclusions. We fit GRB060218, usually considered as an X-Ray Flash (XRF), as a "canonical GRB" within our theoretical model. 
The smallest possible black hole, formed by the gravitational collapse of a neutron star in a binary system, is consistent with the 
especially low energetics of the class of GRBs associated with SNe Ib/c. We give the first evidence for a fragmentation in the fireshell. 
Such a fragmentation is crucial in explaining both the unusually large Tg and the consequently inferred abnormal low value of the 
CBM effective density. 

Key words, gamma rays: bursts — black hole physics — (stars:) binaries: general — stars: neutron 



1. Introduction 
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GRB060218, discovered by the Swift satellite (Masett i et al 
20061) with cosmological re dshift z = 0.033 (Mirab al et al 



2006; Sol lerman et all I2006I) . is the best example of a Gamma- 
Ray Burst (GRB) as sociated with a Supernova (SN) Ib/c 
dCampana et all I2006I) . Its extremely long duration is peculiar, 
with the longest Tgo ever observed (Tgo ~ 2100 s). Tgo is defined 
as the time over which a burst emits from 5% to 95% of its to- 
tal measured energy in the prompt emission. This definition de- 
pends however on the instrumental threshold (see Ruffin fet all 
2007c, for details). This source is also interesting since it repre- 
sents a discrim inant between existing GR B theories: it has been 
pointed out bv lSoderberg et all d2006bl) and lFanetaH d2006l) that 
it is impossible to explain the X- and radio afterglows within 
the traditional synchrotron model. They attempted to fit only the 
late radio data after ~ 10 3 s and they attributed the nature of 
the prompt emission to a yet undetermined "inn er engine" (see 
Soderberg et al., 2006b), possibly a magnetar (Mazzali et al., 



2006). 



In this Letter we present a detailed fit of the entire X- 
and y-ray light curves including the prompt emission: there 
is no need here for the prolonged activity of an inner engine. 
Therefore we explain the unusually high values of the observ ed 
Tgo by our "fireshell" model (see sec. |2j Ruffini et all l2001allR 



I2003L l2005at iBianco & Ruffini 12004 l2005allbl and references 
therein). 

After summarizing our model in sec. |2j in sec. [3] we recall 
GRB060218's observational data. In sec.|4]we show the fit of the 
BAT and XRT light curves (in the 15-150 keV and in the 0.3- 
10.0 keV energy bands respectively, see Figs. 1 1 12b . In Fig.[3]and 
sec. [3] we discuss the actual and effective CircumBurst Medium 
(CBM) density. We outline the occurrence of a fragmentation 
in the fireshell pointing out some analogies with the ejecta of 
Novae. We then proceed to the general conclusions. 

2. The fireshell model 

We assume that all GRBs, whether "short" or "lo ng", originate 
from t he gravitational collapse to a black hole dRuffini et all 
1200 lbl) . The e ± plasma created in the process of the black hole 
formation expands as a spherically symmetric "fireshell" with a 
constant width in the laboratory frame, i.e. the frame in which 
the black hole is at rest. We have only two free parameters char- 
acterizing the source, namely the total energy E'°J of the e ± 
plasma and its baryo n loading B = Mbc 2 /E'°J, where Mb is the 
total baryons' mass (Ruffini et alll2000l) . They fully determine 
the optically thick acceleration phase of the fireshell, which lasts 
until the trans parency condition is reached and the Proper-GRB 
(P-GRB, see Ruffi ni et all 1200 lbl) is emitted. Then, the after- 
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glow emission starts due to the collision between the remaining 
optically thin fireshell and the CBM, and it clearly depends on 
the parameters describing the effective CBM distribution (see 
below). The luminosity of such an afterglow emi ssion consists 
of a r ising branch, a peak, and a decaying tail (Ruffin i et all 
I2001bl) . 

Therefore, u nlike treatments in the current literature (see 
e.g. iPiranl 12004 and references therein), in our model we de- 
fine a "canonical GRB" light curve with two sharply different 
compo nents: the P-GRB and the afterglow (Ruff ini et all 1200 lbl 
2007a). The ratio between the total energies of these two compo- 
nents and the temporal se paration between the ir peaks are func- 
tions of the B parameter (Ruffini et al., 2001b). The peak of the 
afterglow contribute s to what is usually called the GRB "prompt 
emission" (see e.g. iRuffini et all 1200 1 bl 120061) . 

Another crucial assumption is that the afterglow luminos- 
ity is du e to a thermal emission in the co-moving frame of the 
fireshell (Ruff ini et all 120 04). The "R parameter defines the tem- 
perature T of such a thermal emission: 



dE/dt 
4nr 2 crT 4 



(1) 



where dE/dt is the source luminosity, <x is the Stefan-Boltzmann 
constant, r is the radius of the fireshell, A e jj is its effective emit- 
ting area and A r(J is its total visible area. H and the CBM effec- 
tive density n c b m are the two parameters which fully describe the 
effective CBM distribution ta king into account its filamentary 
structure (Ruffinietall|2005b). Similar considerations in a dif- 
ferent context has been recently presented in lPe'er et all d2007l) . 

The description of the engulfment of the CBM matter by 
the fireshell is a most complex and time consuming procedure. 
In the non-relativistic systems such a description can be made 
at each point. In this ultrarelativistic regime a more global ap- 
proach is needed. The arrival time of each photon at the detector 
depen ds on the entire past history of the fireshell (Ruffini et al., 
2001a). All th e observables dep e nds on the eq uitemporal sur- 
faces (EQTSs, iBianco & Ruffinil 12004 l2005al) which produce 
observable effects in a large part of the future history of the 
fireshell. The CBM structure due to its feedback on the equa- 
tion of motion of the fireshell must be therefore inferred self- 
consistently with this very non-linear evolution of the entire 
fireshell history. These difficulties are further increased by the 
necessity to fit the observed light curves in selected energy bands 
(two in the present c ase of GRB060218 and up to five in the case 
of GRB050315, see lRuffini et"aill2006l) . The fulfillment of these 
constraints represents a severe test not only for the validity of the 
theory but also for the spectral models assumed in the data re- 
duction. 

The fireshell mod el has be e n suc cessfully ap- 
plied to GRB050315 (IRuffini et all l2006h. GRB0312q 3 
ilBernardini et"ail 120051). GRB 980425 (IRuffini et all |2007b), 



GRB030329 dBernardinLetaU |2006|). GRB970228 (Bernardin i 
et al., in preparation), GRB991216 IRuffini et all I2003L 120041) . 
Not all the s e sou rces fulfills the correlation proposed by 
Ama ti et ail d2002l) between the isotropic equivalent energy 
emitted in the prompt emission and the peak energy of the 
corresponding time-integrated spectrum (see Dainotti et al., in 
preparation). 



3. GRB060218-SN2006aj 

GRB060218 triggered the BAT instrument of Swift on 18 
February 2006 at 03:36:02 UT and has a T 90 = (2100 + 100) 
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Fig. 1. GRB060218 prompt emission: a) our theoretical fit (blue 
line) of the BAT observations in the 15-150 keV energy band 
(pink points); b) our theoretical fit (red line) of the XRT observa- 
tions in the 0.3-10 ke V energy band (green points) (Data from: 
ICampana etaill2006l) . 



s dCusumano et all l2006l). Th e XRT instrument dKennea et all 
2006; ICusumano et all |2006) began observatio ns ~ 153 s after 
the B AT trigger and continued for ~ 12.3 days feakamoto et all 
2006). The source i s characterized by a flat y-ray light curve 
and a soft spectrum dBarbier et all 120060 . It has an X-ray light 
curve with a long, slow rise and gradual decline and it is con- 
sidered an X-Ray Flas h (XRF) since its pe ak energy occurs 
at E p = 4.9+|]^ keV dCampana et all 120061) . It has been ob- 
served by the Chandra satellite on February 26.78 and March 
7.55 UT (t 8 8 and 17 .4 days) for 20 and 30 ks respectively 
dSoderberg et all l2006bl). The spectroscopic red shift has been 
found to be z = 0.033 dSollerman et all 120061; iMirabal et all 
2006). The correspondi ng isotropic equivalent energy is Ej SO = 
(1.9 + 0.1) x 10 49 erg dSakamoto et all 120061) which sets this 
GRB as a low luminous one, consis t ent wit h most of the GRBs 
associated with SNe dLiang et all l2006bt ICobb et all [2006; 
iGuetta & Delia Valid 120071) . 

GRB060218 is as sociated wi t h SN2006aj whose expansion 
velocity is v ~ 0.1c dPian etail 120061; iFatkhullin et ail 120061; 
ISoderberg et"ail l2006ab ICobbetall 120061) . The host galaxy 
of SN2006aj is a low luminos ity, metal poor star forming 
dwarf galaxy dFerrero et all 12006) with an irregular morphology 
dWiersema et alll2007l). similar to the ones of other GRBs asso- 
ciated with SNe dModjaz et all 120061; ISollerman et alll2006l) . 

4. The fit of the observed data 

In this section we present the fit of our fireshell model to the 
observed data (see Figs.[T]|2]i. The fit leads to a total energy of 
the e ± plasma E'°J = 2.32 x 10 50 erg, with an initial temperature 
T = 1.86 MeV and a total number of pairs N e ± = 1.79 x 10 55 . 
The second parameter of the theory, B = 1.0 x 10~ 2 , is the 
highest value ever observed and is close to the limit for the sta- 
bility of the adiabatic opticall y thick acceleration phase of the 
fireshell (for further details see Ruffin i et al. , 2000). The Lorentz 
gamma factor obtained solving the fireshell equations of motion 
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dBianco & Ruffini l2005allbh is y = 99.2 at the beginning of the 
afterglow phase at a distance from the progenitor r Q = 7.82xl0 12 
cm. It is much larger than y ~ 5 estimated by iKaneko et al.l 
d2006h and lToma et al.1 d2006l) . 

In Fig. Q] we show the afterglow light curves fitting the 
prompt emission both in the BAT (15-150 keV) and in the 
XRT (0.3-10 keV) energy ranges, as expected in our "canonical 
GRB" scenario (see Dainotti et al., in preparation). Initially the 
two luminosities are comparable to each other, but for a detector 
arrival time r a > 1000 s the XRT curves becomes dominant. The 
displacement between the peaks of these two light curves leads 
to a theoretically estimated spectral lag grea ter than 500 s in per - 
fect agreement with the observations (see Lian g et all [2006a). 
We obtain that the bolometric luminosity in this early part coin- 
cides with the sum of the BAT and XRT light curves (see Fig. [2j 
and the luminosity in the other energy ranges is negligible. 

We recall that at ff ~ 10 4 s there is a sudden enhance- 
ment in the radio luminosity and there is an optical luminosity 
dominated by the SN2006aj emission (se e lCampana et al.li 2006; 
ISoderberg et alll2006bllFan et al.Ll2006l) . Although our analysis 
addresses only the BAT and XRT observations, for r > 10 18 
cm corresponding to t% > 10 4 s the fit of the XRT data im- 
plies two new features: 1) a sudden increase of the H factor from 
H = 1.0 x 10" 11 to 'R = 1.6 x 10" 6 , corresponding to a sig- 
nificantly more homogeneous effective CBM distribution (see 
FigEb); 2) an XRT luminosity much smaller than the bolomet- 
ric one (see Fig. |2). These theoretical predictions may account 
for the energetics of the enhancement of the radio and possibly 
optical and UV luminosities. Therefore, we identify two differ- 
ent regimes in the afterglow, one for f l a < 10 4 s and the other 
for ff > 10 4 s. Nevertheless, there is a unifying feature: the de- 
termined effective CBM density decreases with the distance r 
monotonically and continuously through both these two regimes 
from n c bm = 1 particle/cm 3 at r = r a to n c b m = 10" 6 particle/cm 3 
at r = 6.0 x 10 18 cm: n chm oc r a , with 1.0 ^ a ^ 1.7 (see Fig. 

SO- 

Our assumption of spherical symmetry is supported by the 
observations which s et for GRB060218 an opening beaming an- 
gle larger than ~ 37° dLiang et al.Ll2006btlC ampana et al., 2006; 
ISoderberg et all 120061* iGuetta & Delia ValleL 120071) . 

5. The fireshell fragmentation 

GRB060218 presents different peculiarities: the extremely long 
Tgo, the very low effective CBM density decreasing with the dis- 
tance and the largest possible value of B = 10" 2 . These pec uliar- 
ities appear to be correlated. Following Ru ffini et al.l d2007bl) . we 
propose that in the present case the fireshell is fragmented. This 
implies that the surface of the fireshell does not increase any 
longer as r 2 but as ifi with B < 2. Consequently, the effective 
CBM density n c /, m is linked to the actual one n°fL by: 

n c bm = RsheuK.Z m , with <R sheU = (r* / rf , (2) 

where r* is the starting radius at which the fragmentation occurs 
and a — 2 - B (see Fig. [3^). For r* = r a we have n"? = 1 
particl es/cm 3 , as expected for a "canonical GRB" (Ruffini et al., 
I2007ah and in agreement wi th the apparent absence of a mas- 
sive s tell ar win d in the CBM (Soderberg et al., 2006b; IFan et al.L 
l2006l O l2007h . 

The "R parameter defined in Eq.(Q~|) has to take into account 
both the effect of the fireshell fragmentation ("Rsheii) and of the 
effective CBM porosity CR c bm)' 

K = K she u x Kb,n ■ (3) 




Detector arrival time (t^) (s) 

Fig. 2. GRB060218 complete light curves: our theoretical fit 
(blue line) of the 15-150 keV BAT observations (pink points), 
our theoretical fit (red line) of the 0.3-10 keV XRT observations 
(green points) and the 0.3-10 keV Chandra observations (black 
points) are represented together with our theoreti cally computed 
bolometric luminosity (black line) (Data from: ICampana et all 
120061: ISoderberg et al.Ll2006bh . 




Distance (cm) 



Fig. 3. The CBM distribution parameters: a) the effective CBM 
number density (red line) monotonically decreases with the dis- 
tance r following Eq.® (green line); b) the % parameter vs. dis- 
tance. 

The phenomenon of the dumpiness of the ejecta, whose 
measure is the filling factor, is an aspect well known in astro- 
physics. For example, in the case of Novae the filling factor has 
been measured to be in the range 10 2 — 10" 5 (lEderoclite et all 
2006). Such a filling factor coincides, in our case, with 7? v /, ( .//. 

6. Binaries as progenitors of GRB-SN systems 

The majority of the existing models in the literature appeal to a 
single astrophysical phenom enon to explain both the GRB and 
the SN ("collapsar", see e.g. IWooslev & Bloonj, 120061) . On the 
contrary, a distinguishing feature of our theoretical approach is 
to differentiate between the SN and the GRB process. The GRB 
is assumed to occur during the formation process of a black 
hole. The SN is assumed to lead to the formation of a neutron 
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star (NS) or to a complete disruptive explosion without rem- 
nants and, in no way, to the formation of a black hole. In the 
case of SN 2006aj the formati on of such a NS has been actually 
inferred by IMaedae t al. (2007) because of the large amount of 
58 Ni (0.05 M ). Moreover the significantly small initial mass of 
the SN progenitor star M « 2OM is expected t o form a NS 
rather than a black hole when its core collapse s dMaeda et all 
20071: iFerrero et all 120061: iMazzali et all 120061: iNomoto et ail 
2007). In order to fulfill both the above requirement, we as- 



sume that the progenitor of the GRB and the SN consists of 
a binary system formed by a NS close to its critical mass col- 
lapsing to a black hole, and a companion star evolved out of the 
main sequence originating the SN. The temporal coincidence be- 
tween the GRB and the SN phenomenon is expl ained in term 
of the concept of "induced" gravitational collapse (Ruffin iet all 
2001c H2007bl) . There is also the distinct pos sibility of observing 
the young born NS out of the SN (see e.g. IRuffini et all l2007bl 
and references therein). 

It has been often proposed that GRBs associat ed with SNe 
Ib/c, a t smaller redshift 0.0085 <z< 0.168 (see e.g. lDella VaiR 
2006, and references therein), form a different class, less lumi- 
nous and possibly much mor e numerou s than t he high luminos- 
ity GRBs at higher redshi f t dPian et all 120061: ISoderberg et all 
2004; IMaedaet all l2007t iDella ValleL 120061) . Therefore they 
have been proposed to originate from a separ ate class of progen- 
itors dLiang et all l2006bt ICobb et all 120061) . In our model this 
is explained by the nature of the progenitor system leading to 
the formation of the black hole with the smallest possible mass: 
the one formed by the colla pse of a just overcritical NS (Ruffini, 
120061: IRuffini et all l2007bl) . 

The recent ob servation of G RB06061 4 at z = 0.125 withou t 
an associated SN (Delia Valle et al l l2006l:lMangano et alll2007h 
gives strong support to our scenario, alternative to the collapsar 
model. Also in this case the progenitor of the GRB appears to 
be a binary system composed of two NSs or a NS and a white 
dwarf (Caito et al., in preparation). 

7. Conclusions 

GRB060218 presents a variety of peculiarities, including its ex- 
tremely large Tgo and its classification as an XRF. Nevertheless, 
a crucial point of our analysis is that we have successfully ap- 
plied to this source our "canonical GRB" scenario. 

Within our model there is no need for inserting GRB060218 
in a new class of GRBs, such as the XRFs, alternative to the 
"canonical" ones. This same point recently receive d strong ob- 
servat ional support in the case of GRB060218 dLiang et all 
2006 a]) and a cons ensus by other models in the literature 
( iKanekoet ail 120061) . 

The anomalously long Tgo led us to infer a monotonic de- 
crease in the CBM effective density giving the first clear ev- 
idence for a fragmentation in the fireshell. This phenomenon 
appears to be essential in understan ding the features of also 
other GRBs (see e.g. GRB050315 in IRuffini et all l2007bl and 
GRB970228 in Bernardini et al., in preparation). 

Our "canonical GRB" scenario originates from the gravita- 
tional collapse to a black hole and is now confi rmed over a 10 6 
range in energy (see e.g. IRuffini et all l2007al and references 
therein). It is clear that, although the process of gravitational col- 
lapse is unique, there is a large variety of progenitors which may 
lead to the formation of black holes, each one with precise signa- 
tures in the energetics. The low energetics of the class of GRBs 
associated with SNe, and the necessity of the occurrence of the 
SN, naturally leads in our model to identify their progenitors 



with the formation of the smallest possible black hole originat- 
ing from a NS overcoming his critical mass in a binary system. 
For GRB060218 there is no need within our model for a new or 
unidentified source such as a magnetar or a collapsar. 

GRB060218 is the first GRB associated with SN with com- 
plete coverage of data from the onset all the way up to ~ 10 6 
s. This fact offers an unprecedented opportunity to verify theo- 
retical mo dels on such a GRB class. For example, GRB060218 
fulfills the lAmati et all (2002) relation unlike other sources in its 
same class. This is particularly significant, since GRB060218 is 
the only source in such a class to have an excellent data coverage 
without gaps. We are currently examining if the missing data in 
the other sources of such a class may have a prominent role in 



their non-fulfillment of the Amati et al. ( 20021) relati on (Dainotti 
et al., in preparation; see also Ghisellini et all [20061) . 
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